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Abstract 

Aims. The purpose of this study is to investigate the dis- 
tributions of the isomeric molecules HCN and HNC and 
estimate their abundance ratio in the protostellar core Cha- 
MMS 1 located in Chamaeleon i. 

Methods. The core was mapped in the 7=1-0 rotational 
hnes of HCN, HNC, and HN'^C. The column densities of 
H'^CN, HN'^C, H'^NC and NH3 were estimated towards 
the centre of the core. 

Results. The core is well delineated in all three maps. 
The kinetic temperature in the core, derived from the 
NH3 (1,1) and (2,2) inversion Hnes, is 12.1 + 0.1 K. The 
HN'^C/H'^CN column density ratio is between 3 and 4, i.e. 
similar to values found in several other cold cores. The 
HN'^C/H'^'NC column density ratio is ~ 7. In case no 
'""N fractionation occurs in HNC (as suggested by recent 
modelling results), the HNC/HN'^C abundance ratio is in 
the range 30 - 40, which indicates a high degree of '^C 
fractionation in HNC. Assuming no differential '^C frac- 
tionation the HCN and HNC abundances are estimated to 
be ~ 7 10"'" and ~ 2 10"^, respectively, the former being 
nearly two orders of magnitude smaller than that of NH3. 
Using also previously determined column densities in Cha- 
MMSl, we can put the most commonly observed nitroge- 
nous molecules in the following order according to their 
fractional abundances: ;t'(NH3) > ;i'(HC3N) > ;i'(HNC) > 
;r(HCN)>;^(N2H+). 

Conclusions. The relationships between molecular abun- 
dances suggest that Cha-MMSl represents an evolved 
chemical stage, experiencing at present the Tate-time' 
cyanopolyyne peak. The possibility that the relatively high 
HNC/HCN ratio derived here is only valid for the '^C iso- 
topic substitutes cannot be excluded on the basis of the 
present and other available data. 



'Based on observations collected at the European Southern Observa- 
tory, La Silla, Chile, and at the Parkes Observatory operated by the Aus- 
tralia Telescope National Facility (ATNF). 



1 Introduction 

Hydrogen cyanide, HCN, and its metastable isomer 
hydrogen isocyanide, HNC, are commonly used trac- 
ers of dense gas in molecular clouds. The HNC 
column density has been observed to be similar or 
larger than that of HCN in cold regions, whereas in 
warm GMC cores the [HNC]/[HCN] ratio is typically 
much smaller than unity (e.g., I Goldsm ith et al. 19 811: 

fChurchwell et al. 1984|; ^rvine & S chloerb 19841 : 

fSchilke et al. 1992J : LHirotaetal. 199 81). 

Gas-phase chemistry models, including both ion- 
molecule and neutral-neutral reactions, predict that the 
HNC/HCN abundance ratio should be close to unity in 
cold gas (see e.g. |Herbstetal. 2000), where also pre- 
vious work is reviewed). The two isomers are thought 
to form primarily via the dissociative recombination of 
HCNH^, which in turn is produced by a reaction be- 
tween NH3 and C^. The recombination reaction regulates 
the abundance ratio despite the initial production mech- 
anisms of the two isomers, provided that they are effi- 
ciently protonated via reactions with ions, such as H^ and 
HCO+ (I Churchwelletal. 19841 : IBrown et al. 19 891). Re- 
cent theoretical studies corroborate the equal branching 
ratios for HCN and HNC in dissociative recombination 
( IHickman et al. 20051 : llshii et al. 2006 11 

HNC/HCN column density ratios clearly larger than 
unity have been reported: IChurchwelletal. 19841 find 
that [HNC]/[HCN] may be as high as ~ 10 in some dark 
clouds cores, and |Hirota et al. 1998 1 find a ratio of ~ 5 in 
L1498. The first result quoted is uncertain because it de- 
pends on the assumed degree of '^C fractionation in HNC, 
whereas the second value refers to HN'^C and H'^CN 
measurements. Nevertheless, both results suggest that the 
HCN/HNC chemistry is not fully understood. Further deter- 
minations of the relative abundances of these molecules to- 
wards cores with well-defined chemical and physical char- 
acteristics seem therefore warranted. 

In this paper we report on observations of HCN and 
HNC and some of their isotopologues towards the dense 
core Cha-MMSl, which is situated near the reflection neb- 
ula Ced 110 in the central part of the Chamaeleon 1 dark 
cloud. We also derive the gas kinetic temperature in the core 
by using the (1,1) and (2,2) inversion lines of NH3. Cha- 
MMSl was discovered in the 1.3 mm dust continuum by 
I Reip urth et al. 199'6|, an d was studied in several molecules 
by I Kontinen et al. 20001 . The object is embodied in one 
of the most massive 'clumps' in Chamaeleon i identi- 
fied in the large scale C'^O survey of IHaikala et al. 20051 
(clump no. 3, ~ 12Mo). |Reipurth et al. 1996| sug- 
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gested that Cha-MMSl contains a Class protostar. A 
FIR source (Ced 110 IRS 10) was detected near the cen- 
tre of the core by ILe htinen et al. 2001 J . However, no 
centimetre continuum nor near-infrared sources have been 
found in its neighbourhood (TLehtinen et al. 200 31 and ref- 
erences therein), and the core therefore represents a very 
early stage of star formation. At the same time, its chem- 
ical composition has probably reached an advanced stage 
r iKontinen et al. 200 01). 

In Sect. 2 of this paper we describe the observations. In 
Sect. 3 we present the direct observational results in the 
form of spectra and maps, and derive the H'^CN, HN'^'C, 
and H'^NC column densities. By combining our results 
with previous observations we derive in Sect. 4 the frac- 
tional abundances of several nitrogen containing molecules 
in Cha-MMSl, and discuss briefly the chemical state of the 
core. In Sect. 5 we summarize our results. 



2 Observations 



molecule ° 

(Debye) 



transition 



(MHz) 



The 7=1-0 transitions of HCN, H'^CN, HNC, HN'"C and 
H'^NC at about 90 GHz were observed with the Swedish- 
ESO-Submillimeter- Telescope, SEST, located on La Silla 
in Chile. The observations took place in December 1990. 
A 3 mm dual polarization Schottky receiver was used in 
the frequency switching mode. The backend was a 2000 
channel acousto optical spectrometer with an 86 MHz band- 
width. The velocity resolution with this configuration is 
about O.lSkms ' at 90 GHz. Typical values for the sin- 
gle sideband system temperatures ranged from 350 K to 
450 K. The half-power beam width of the SEST is about 
55" at the frequencies observed. Calibration was achieved 
by the chopper wheel method. Pointing was checked every 
2-3 hours towards a nearby SiO maser source. We estimate 
the pointing accuracy to have been better than 5" during the 
observations. The focusing was done using a strong SiO 
maser 

The (1,1) and (2,2) inversion lines of NH3 at about 23 
GHz were observed with the Parkes 64-m radio telescope 
of the Australian Telescope National Facility (ATNF) in 
April 2000. The FWHM of the telescope is about 80" at 
this frequency. The two lines were observed simultaneously 
in two orthogonal polarizations in the frequency switching 
mode, using the 8192 channel autocorrelator The calibra- 
tion was checked by observing standard extragalactic cal- 
ibration sources at different elevations, and by comparing 
the ammonia line intensities towards some strong galactic 
sources also visible from the Eifelsberg 100-m telescope, 
taking the diff'erent beam sizes into account. 

The electric dipole moments of the molecules and the line 
frequencies used in this paper are listed in Tabled Except 
for the frequencies of the HNC and HN'^C hyperfine tran- 
sitions adopted from |Frerking et al. 1979| , the parameters 
are obtained from the Jet Propulsion Laboratory molecular 
spectroscopy data base (http://spec.jpl.nasa.gov), 
where also the original references are given. The hyper- 
fine components of NH3 are not listed. The frequencies 
given represent the centre frequencies of 18 and 21 indi- 
vidual components of the (1,1) and (2,2) lines, respectively, 
distributed over a range of about 4 MHz. 
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* Frequencies from |Frerking et al. 1979) 

Table 1 : The permanent dipole moments of the molecules 
and the rest frequencies of the observed transitions. 

3 Maps and column densities 

3.1 Maps 

The integrated HCN(7 = 1-0), HNC(7 =1-0) and 
HN'^C(7 =1-0) intensity maps of Cha-MMSl are pre- 
sented as contour plots in Figs.m2land|3l respectively. The 
values of the contour levels are given in the bottom right of 
each figure. The position of the 1.3 mm dust emission max- 
imum, Cha-MMSla (| Reipurth et al. 1996J ), is chosen as 
the map centre. The coordinates of Cha-MMSl a are aaooo = 
11*'06'"31.7\ ^2000 = -77°23'32". The grid spacing used 
in the maps was 30". Also indicated in these figures are 
the location of the HN'^C maximum (Aa, A6 = -20", 20", 
black dot), and the positions used for the HCN(7 =1-0) 
excitation temperature estimates in the LTE method (open 
circles, see below). Long integrations towards the HN'-'C 
maximum were made besides HNC and HCN also in the 
7=1-0 lines of H'^CN and H'^NC. The five spectra 
obtained towards this position are shown in Fig.|5l 

The three maps have differences, but they all show a rel- 
atively compact distribution peaking near the position of 
Cha-MMSl. The HN^^C emission with the lowest opacity 
reflects probably best the distribution dense gas. 

3.2 LTE column densities 

The column densities of the molecules have been deter- 
mined with two different methods. First, we have derived 
the optical thicknesses of the lines and the total column den- 
sities of the molecules by assuming local thermodynamic 
equilibrium, LTE. The LTE method involves idealistic as- 
sumptions about the homogenity of the cloud. In the next 
subsection we present modelling of the observed spectra 
with a Monte Carlo program using a realistic density dis- 
tribution. 

In the following, we assume that the rotational levels of 
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Figure 1: Integrated HCN(7 - 1 -0) line intensity map. The 
integration is performed over all three hyperfine structure 
components {F-l-l,F = 2- l and F - 0-1). The 
map is centered on the position of Cha-MMSla (denoted 
by a cross), ajooo = 11^06™31.7\ 52om = -77°23'32", and 
the offsets are in arcseconds. The contour levels indicated 
with ten shades of grey and their values can be read in the 
bottom right corner. The intensity unit is K km s The 
positions used for the column density estimates (-20", 20") 
and Tex estimates (see TableO are marked with a black dot 
and open circles, respectively. The grid spacing used in the 
map is 30". 
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Figure 3: Integrated HN'^C(7 =1-0) line intensity map. 
The coordinates and markers as in Fig^ 

Aa Ad Tex Tmain Au A^(HCN) 

(") (") (K) (kms-') (cm'^) 

20 1.1 +0.3 0.49 + 0.30 4.7 10'^ 

-50 -10 5.3!^-(^, 0.9 +0.4 0.58 + 0.07 4.6 lO'^ 

10 -10 5.6!f° 0.7 +0.3 0.56 + 0.08 3.8 lO'^ 

70 20 5.4;j;° 0.5 +0.3 0.49 + 0.04 2.010'^ 

Table 2: Estimates of the HCN(7 =1-0) excitation tem- 
peratures and the HCN column densities towards selected 
positions using hyperfine structure fits. The coordinates are 
offsets from the position of Cha-MMS 1 a. 
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Figure 2: Integrated HNC(y =1-0) line intensity map. 
The coordinates and markers as in Fig [2 



all five isotopologues are characterized by the same Tex, and 
that the core is homogeneous, i.e. the same Tex is valid for 
different locations. The former assumption can be justified 
by the similar dipole moments and rotational constants of 
the molecules. The second assumption infers that the exci- 
tation conditions remain constant along the line of sight and 
also between close-lying positions in the core. 

The excitation temperature is derived by fitting first the 
optical thickness, t, to the observed hyperfine ratios of a 
"clean" HCN spectrum, i.e. a spectrum with a good S/N 
and hyperfine intensity ratios in accordance with the LTE 
assumption. The Tex is solved from the antenna equation us- 
ing observed peak main beam brightness temperature, Tmb 
and assuming uniform beam filling. The fit is possible to- 
wards the four positions listed in Table|2 In other positions 
with strong signal, like the HN'-'C maximum (-20", +20"), 
or the HCN maximum at (-20", -20") the HCN spectra 
show 'anomalous' hyperfine ratios. The smallest error is 
obtained towards the offset (40", 20") near the edge of the 
core with Tex ~ 5 K. This HCN spectrum is shown in Fig.|3] 
The other positions listed in Table|2]give similar values of 
Tex but the errors are larger. Also HCN column density esti- 
mates using the values for r and At; obtained from hyperfine 
fits are given in Table |2l 

The HCN, H^^CN, HN^^C and H'^NC column densities 
towards the HN'^C maximum, (-20", +20"), have been de- 
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Rest Frequency (MHz) 
8.8638 10 8-8636 10 




Figure 4: The HCN(7 =1-0) spectrum towards offset 
(40", 20") used to determine the T^^. The dotted line repre- 
sents a fit to the hyperfine structure. 
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Table 3: Column densities towards the HN'-'C maximum 
derived using the LTE method. The assumed value for T^x is 
5.2;^Q g K. ° The intensity value is for the weakest hyperfine 
component F -Q)-\ which comprises one ninth of the total 
intensity of the 7=1-0 line. A column density estimate 
is not given because of the obvious large optical thickness. 
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Figure 5: Observations used for the column density esti- 
mates. The position corresponds to the HN^"^C maximum, 
which lies at the offset (-20", +20") from the dust contin- 
uum peak Cha-MMSla (marked with a black dot in Figs.Q} 



rived from the integrated 7=1-0 line intensities and 
the Tex estimate derived from the HCN spectrum shown in 
Fig. 13 For HCN the weakest of the well separated hyperfine 
components, F = - 1, with the least optical thickness has 
been used in this calculation by taking its relative strength 
into account. For HNC we did not attemp to do a simi- 
lar estimate because the hyperfine components overlap and 
are very likely optically thick. Overlapping hyperfine com- 
ponents hamper also the HN'-'C column density estimate, 
but their optical thicknesses should be lower than those of 
HNC by an order of magnitude or more. H'^NC has no hy- 
perfine structure due to zero nuclear spin of '^N. The results 
of these calculations are presented in Table|3 

The following column density ratios can be derived from 
Tabled 



even more obvious in the case of the HNC/HN^C intensity 
ratio, 1.7. We do not use therefore the values derived from 
the HCN and HNC spectra towards the core centre. The 
HCN and HNC spectra are discussed at the end of Sect. 4. 

The ammonia column density, A^(NH3), and the gas ki- 
netic temperature, r^in, were derived using hyperfine fits to 
the NH3 inversion lines arising from the rotational states 
{J,K) - (1,1) and (2,2). The spectra are shown in 
Fig. |6l The standard method of analysis described e.g. 
in I H o et al. 19791 gave the following values: A^(NH3) = 
(1.4 + 0.2) 10'^ cm-2 and Tun = 12.1 ± 0.1 K. Here it is 
assumed that the relative populations of all metastable lev- 
els of both ortho- and para-NHa are determined by thermal 
equilibrium at the (2, 2)/(l, 1) rotational temperature. 



A?[HNi3C] 
A^[Hi3CN] 
A?[HNi3C] 
A?[Hi5NC] 
M[HCN] 
Ar[Hi3CN] 



3.3 +0.4 

7.4 + 1.8 
17 + 2 . 



The HCN/H'^CN column density ratio, 17, is much 
lower than the isotopic '^C/'^C ratios usually observed in 
nearby molecular clouds or local diffuse ISM (~ 60, e.g. 
(Lucas & Liszt 1998J; the terrestrial value is 89), and is 
likely to be affected by absorption or scattering in the less 
dense envelope of the core. Radiative transfer effects are 



3.3 Monte Carlo modelling 

H'^NC, HN'^C, HCN and H'^CN spectra were calcu- 
lated with the Monte Carlo radiative transfer program 
of lJuvela 19971 . A constant kinetic temperature and 
a spherically symmetric density distribution, n~ r 
were assumed. The selected density power law conforms 
with observational results towards star-forming cores (e.g., 
I Tate matsu et al. 20041 ). In order to avoid divergence in the 
cloud centre, we assume a constant density within a dis- 
tance corresponding to 6% of the outer radius of the model 
cloud. Compared with the beam size, the constant density 
region is small and the results are not very sensitive to the 
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Figure 6: NUi(J,K) 
towards Cha-MMSla. 
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(1,1) and (2,2) inversion spectra 



selected value of its radius. The spectra were convolved 
with a gaussian that corresponds to the size of the beam 
used in the observations. The simulated spectra were fitted 
to observations by adjusting four parameters: the outer ra- 
dius of the core, the scaling of gas density, and the frac- 
tional abundances of H'^CN and HN'^C. The HCN fit was 
optimized for the weakest hyperfine component F = - 1. 
The abundances of the other species were determined by as- 
sumed fixed ratios ['^C]/['-'C]=20, as suggested by the LTE 
column density determination, and [''*N]/['^N]=280 (terres- 
trial). The resulting fractional H'^CN and HN'^C abun- 
dances do not, however, depend on these assumptions about 
isotopic ratios. The calculations were done with rkin=12K, 
which was derived from the ammonia observations, and 
=20 K.The lower temperature resulted in better fit although 
the difference in the;^'^ values was no more than 25%. The 
results for the 12 K are shown in Fig.0 The outer radius 
of the model is 60". The volume density increases from 
2.3 10"* cm"-' at the outer edge to 1.4 10^ cm"-' in the cen- 
tre, and the molecular hydrogen column density towards the 
centre is A?(H2) = 1 .9 lO^^ ^,^-2 

The fractional HN'-'C and H'-'CN abundances resulting 
from the fit are 2.1 10"'" and 5.5 10"", respectively. The 
HN'-'C/H'-'CN abundance ratio obtained from the simula- 
tions is thus 3.8, i.e. close to the value obtained from the 
LTE method. 

The largest discrepancy between the observed and mod- 
elled profiles is seen in HCN. It can be explained if 
the optical depth is larger than in this particular model. 
For example, a foreground layer of cold gas would read- 
ily decrease the intensity of the HCN lines and would 
bring the ratios between the hyperfine components closer 
to the observed values. The hyperfine ratios are, how- 
ever, sensitive to the actual structure of the source 
(I Gonzales- Alfonso & CernicharoJ993 1). 

The optical thickness of the HN'^C main component is 
~ 0.7. This moderate value suggests that the method used 
in the previous subsection results in a slight underestimata- 
tion of column densities. The H'^CN and HN'^C column 
densities (fractional abundances multiplied by A^(H2)) are 
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Figure 7: Simulated spectra fitted to the observations to- 
wards the HN"C maximum. The assumed kinetic temper- 
ature, 7'kin=12K, is based on ammonia observations. The 
following isotopic ratios were used in the model calcula- 
tions; [i2c]/[i'C]=20, and [''*N]/[15N]=280. 



indeed about two times larger than those presented in the 
previous section: A?(Hi3CN) = 1.1 lOi^cm"^ A^(HN"C) = 
4.0 lO'^cm"^. These results suggest, however, that the LTE 
method does not lead to crude underestimates of the column 
densities of rarer isotopologues of HCN and HNC. 



4 Discussion 

4.1 Fractional abundances 

The present observations and the LTE method give reason- 
able estimates for the column densities of H'^CN, HN'^C, 
H'^NC and NH3 neai- the centre of Cha-MMS 1 . Unfortu- 
nately the column densities of the common isotopologues 
of HCN and HNC cannot be directly determined because 
of the large optical thickness of their lines. H'^NC is, 
however, likely to be useful for this purpose. Accord- 
ing to the modelling results of ITerzieva & Herbst 20001 
hardly any '^N fractionation occurs in dense interstellar 
clouds. The isotopic '''N/'^N ratio lies probably between 
240 and 280, i.e. the values derived for the local ISM 
CILuc as & Liszt 199 81). and the terrestrial atmosphere. Us- 
ing this range and the H'^NC column density listed in Ta- 
ble 13 we obtain A^(HNC) ~ 610'^cm"2. Assuming that 
no differential '^C fractionation occurs, we get from the 
HN'^C/H'^CN ratio derived above that A?(HCN) ~ 2 10"^ 



-2 

cm . 

The dust continuum emission provides the most reliable 
estimate for the H2 column density. [Reipurth et al. 1996 J 
measured towards Cha-MMS 1 a total 1.3 mm flux den- 
sity of 950 mJy from a nearly circular region with a di- 
ameter of about 50" (bordered by the lowest white con- 
tour in Fig. 8). As this size is similar to the SEST beam- 
size at 3 mm, the resulting average intensity and average 
A^(H2) can be compared with the column densities derived 
here. The dust temperature as derived from ISO observa- 
tions by ILehtinen et al. 20011 is 20 K (Cha-MMS 1 = Ced 
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1 10 IRS 10). Due to crowding of sources in this region and 
the Hmited angular resolution of ISOPHOT, this value may 
be affected by emission from other sources in its vicinity 
(see Fig. 2 of IILehtinen et al. 200n '). Therefore, value 20 
K is probably an upper limit. A lower limit is provided by 
the gas kinetic temperature ~ 12 K derived above. By us- 
ing the two temperatures and Eq. (1) of IMotte et al. 19981 
with /fijmm - 0.01 gcm"^ appropriate for circumstellar 
envelopes around Class and Class 1 protostars, we get 
A^(H2) - 1.9 - 4.010^^cm"^. The corresponding mass 
range of the dust core obtained from the total flux given 
by |Reipurtheta~1996J is 0.4 - 0.8 Mq. 

The median of the reasonable H2 column density 
range, A^(H2) - S.OlO^^cm"^, yields the following val- 
ues for the fractional HCN, HNC and NH3 abundances: 
;if(HCN) ~ 7 10-"', Ar(HNC) ~ 2 10"^ and xi^^i) ~ 
5 10"^. The column densities of two other nitrogen bear- 
ing molecules towards Cha-MMSl were determined by 
tKontine n~al. 20001 : A^CHCsN) = (4.5 ± 0.8) lO'"* cm^^, 
and A^(N2H+) = (1.4 + 0.1) 10'^ cm^^, which imply the 
fractional abundances ;if(HC3N) ~ 210"^ and;i'(N2H+) ~ 
5 10-1°. 

We find that both HCN and N2H+ are about 100 times 
less abundant than NH3, whereas HC3N is almost equally 
abundant as NH3. The derived N2H^/NH3 column den- 
sity ratio is in accordance with previous observational re- 
sults (e.g., ITa fallaetal. 2 0021; IHotzel et al. 20041 ). The 
abundance ratios are discussed in Sect. 4.3. The C"^0 
column density in this direction is (1.5 + 0.1) 10'^ cm"^ 
( IHaikalaet al. 20051 ). The fractional C"^0 abundance is 
therefore 5 10"^ or smaller This upper limit falls below the 
'standard' value ~ 2 10"^ of [ IFrerking et al. 1982| , and sug- 
gests CO depletion. 

4.2 Comparison with previous maps 

The superposition of our HN'^C map, the HC3N map of 
(Kontinen et al. 2000 1, and the 1.3 mm dust continuum map 
of |,Reipurth et al. 1996 J is shown in Fig.lH The two molec- 
ular line maps have similar angular resolution (about 60") 
whereas the resolution of the continuum map is 22". Both 
molecules show rather compact, symmetric distributions 
around the dust continuum source, HC3N agreeing slightly 
better with dust than HN'^C. Nevertheless, the dense core 
can be clearly localized with the aid of HN^-^C as well as 
HCN and HNC. This is in contrast with SO, which shows 
two peaks on both sides of the core ( | Kontin en et al. 20 00 1 ), 
and with C"^0 which has a flat distribution in this region 
dHaikal a et al. 20051 ). probably due to depletion. This 
result agrees with the results of IHirota et al. 2003*1 and 
INikolic etal. 20031 . which suggest fliat HCN and HNC are 
less affected by the accretion onto dust grains than some 
other commonly observed molecules. 

4.3 The chemical state of Cha-MMSl and the 
HNC/HCN abundance ratio 

CO depletion and a large NH3 abundance point towards 
an advanced stage of chemical evolution. The NH3/N2H^ 
abundance ratio ~ 100 agrees with model predictions for 
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Figure 8: The HN C map (dark contours) superposed on 
the HC3N map from IK ontinen et a l. 20001 (grey scale) and 
the 1.3 mm dust continuum map of [Reipurth et al. 1996J 
(white contours). 

dense cores, but does not put severe constraints on the 
evolutionary stage ( lAikawa et al. 20051 ). According to 
the observational results of IHotzel et al. 20041 this value 
should be more appropriate for a star-forming than for a 
prestellar core. However, due to slightly different beam- 
sizes the accuracy of the present determination of the 
NH3/N2H^ ratio is not sufficient for distinguishing between 
the two cases. | Kontinen et al. 2000 1 suggested that the 
large HC3N abundance observed in Cha-MMSl indicates 
either a so called 'late-time' cyanopolyyne peak associated 
with a large degree of depletion according to the model of 
I Ruf fle et al. 19971 . or carbon chemistry revival by the in- 
fluence of a newly born, so far undetected star A large am- 
monia abundance can be understood in either case. Ammo- 
nia production benefits from the diminishing of CO and can 
resist depletion longer (e.g. fNejad & Wagenblast 19991). 
On the other hand, ammonia can possibly form on grain 
surfaces and be released into the gas phase in shocks 
(INej adetal. 19 901; FAikaw a et al. 20031 ). As there is no 
clear evidence for star-cloud interaction in Cha-MMSl, the 
alternative that the core - or the portion traced by the present 
observations - represents matured pre-stellar chemistry is 
more likely. 

The behaviour of HCN and HNC under these circum- 
stances is not clear The main gas-phase production path- 
way via the HCNH^, which is effective especially at later 
stages, requires NH3 and C^ ions. The latter can be pro- 
duced also in deep interiors of molecular clouds via the cos- 
mic ray ionization of C atoms, or in the reaction between 
CO and the He^ ion, the latter mechanism becoming less 
important with CO depletion (|Nejad & Wagenbl ast 1999| ). 
To our knowledge the production of HCN and HNC on in- 
terstellar grains and their desorption has not been studied. 

The models of fAikawa et al. 20051 for collapsing 
prestellar cores include the gas-phase production of HCN, 
HC3N, NH3, and N2H+. Although the order of fractional 
abundances, a'(NH3) > ;(f(HCN) > xi^i'H^), and the 
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NHs/NaH^ abundance ratio derived from the present obser- 
vations are well reproduced by this model, especially for the 
gravity-to-pressure ratio a - 4 impying a rapid collapse and 
less marked depletion, the predicted HCN column densities 
relative to NH3 and N2H^ are an order of magnitude larger 
than those observed here. A likely explanation is that HCN 
is not as centrally peaked than NH3 and NaH^ (see Fig. 3 
of lAikawa et al. 20 051'). Therefore a large fraction of HCN 
resides in on the outskirts of the core, where the density is 
not sufficient for the collisional excitation of its rotational 
levels, and the molecule is not seen in emission (see also 
the next section). Another discrepancy is that the model 
predicts much lower HC3N column densities relative to am- 
monia than observed here. The agreement is slightly better 
for the model with overwhelmingly large gravitational po- 
tential (a - 4) than for the slowly evolving model near hy- 
drostatic equilibrium (a = 1.1) which results in very high 
degrees of depletion. Unfortunately the spectral resolution 
and S/N ratio of the present data do not allow to study the 
core kinematics. 

As compared to the 'principal' N-bearing molecules NH3 
and HC3N, HCN and HNC have clearly smaller abun- 
dances, and probably only marginal importance to the ni- 
trogen chemistry in the core. Nevertheless, the fact that 
HNC/HCN abundance ratio is once again observed to be 
larger than unity remains a problem. ITalbi & Herbst 19981 
have shown that even though the reaction between NH3 
and C^ may produce a significant amount of the metastable 
ion H2NC^, which should produce only HNC in electron 
recombination, the energy released in the former reaction 
leads to the efficient transformation into the linear isomer 
HCNH^ before the recombination. Likewise, according to 
fHerbst et al. 2000.1 . the vibrational energies of HNC and 
HCN formed in the neutral reactions C + NH2 —> HNC -1- H, 
and N + CH2 — > HCN + H are sufficient to overcome the 
isomerization barrier, and roughly equal amounts of both 
isomers are produced. 

The theoretical results presented in 
ITalbi & Herbst 1998J and tHerbsteta l. 20001. and more 
recently by IHickman et al. 201751 and |Ishiietal. 2006 1, 
show quite convincingly that the HNC/HCN abundance 
ratio resulting from the chemical reactions known to be 
effective in dark clouds cannot differ significantly from 
unity. There are two possible explanations for the discrep- 
ancy between the observations and theoretical predictions. 
Either so far unknown processes favour HNC in the cost of 
HCN, or then the derived HN'^C/H'^'CN column density 
ratio is affected by different degrees of '^C fractionation in 
HNC and HCN. In the latter case '^C would about three 
times more enhanced in HNC than in HCN in Cha-MMS 1 . 

The HN'^C/H'^NC column density ratio 7.4 together 
with the above quoted range of '"^N/'^N ratios imply that the 
HNC/HN'^C ratio in Cha-MMSl is in the range 30-40. As 
the typical '^C/'^C isotopic ratio observed in nearby clouds 
is 60 ( ILucas & Li szt 1998 1) the result suggests consider- 
able '^C enhancement in HNC. An HC'^N column density 
determination would readily show, according to the results 
of [Terzieva & Herbst 2000 1, whether also HCN has a sim- 
ilar degree of '^C fractionation. 

In the only survey so far including all these iso- 



topologues, llrvine & Schloerb 19841 found very similar 
HN'^C/H'^NC and H'^CN/HC'^N intensity ratios (~ 1.5) 
towards TMC-1. This observation suggests similar '■'C 
fractionation in HCN and HNC. On the other hand, this 
result might not be univer sally vahd. Ilkeda e t al. 20021 
determined substantially different H'^'CN/HC'^N col- 
umn density ratios in the starless dark cloud cores 
L1521E, TMC-1, and L1498, listed in order of increasing 
H'^CN/HC'^N ratio (2.6, 6.7, and 11). The HN^^C/Hi^CN 
column density ratios in these sources increase in the same 
order (0.54, 1.7, and 4.5; FHirot a et al. 199 81). The val- 
ues obtained towards TMC-1 are consistent with those of 
llrvin e & Schloerb 1984ll . but the results towards L1521E 
and L1498 suggest (if the possibility of '^N fractionation 
is neglected) that the '^C fractionation in HCN changes 
from source to source, and perhaps the variation in the 
case of HNC is even larger L1521E is suggested to rep- 
resent youthful chemistry with little molecular depletion 
( IHirotaetal. 20021 : |Tafalla & Santiago 2004| ) whereas 
the position observed in L1498 is known to be heavily 
depleted ([ Willacy et al. 1998J ). Thus, it appears that the 
freeze-out of molecules favour the '^C substituted isotopo- 
logues of HCN and HNC. It seems warranted to further 
examine the possibility of their differential '^C fractiona- 
tion, although no theoretical arguments are in favour of dif- 
ferent gas-phase production and destruction rates of HN'^C 
and H'^CN, or their different freezing rates. As no empir- 
ical data on this matter is currently available, the best way 
to proceed is probably to perform similar observations as 
those done by llrvine & Schloerb 1984 1 towards a sample 
of cores with different degree of depletion. 

4.4 Scattering envelope 

The low HCN/H'^CN and HNC/HN^^C intensity ratios and 
the line shapes in the HCN and HNC spectra observed to- 
wards the centre of the core require some attention. The 
intensity ratios suggest very large optical thicknesses for 
the HCN and HNC lines, and consequently, deep self- 
absorption features and 'anomalous' hyperfine component 
ratios are to be expected. A clear absorption feature can be 
seen in the simulated HCN spectrum shown in Fig. al- 
though a very low '^C/'^C isotopic ratio is used in order to 
avoid inverted hyperfine ratios, which would disagree with 
the observations. The observed HCN and HNC spectra are, 
however, free from absorption dips, and they look more like 
weakened 'normal' lines. 

This weakening can probably be attributed to less dense 
gas surrounding the core. The HCN and HNC molecules 
in the envelope, where the density lies below the criti- 
cal density of their 7=1-0 lines, function as scat- 
terers of line emission arising from the core, i.e. they 
re-emit absorbed photons (into random directions) before 
colliding with other particles. The effect is less marked 
for '^^C isotopologues, and the resulting HCN/H'^CN and 
HNC/HN'^C intensity ratios decrease, without HCN and 
HNC loosing their characteristic line patterns. This model 
was used by ICernicharo etaL 1 9841 to explain HCN hy- 
perfine anomalies, particularly the intesity of the intrinsi- 
cally weakest component F = - 1 observed towards dark 
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clouds. The HCN spectrum shown in Fig.|5]conforms with 
this explanation: the F - - 1 component on the left is 
stronger than F - 1 - I component on the right, although 
the optically thin LTE ratio is 1 : 3. 

The scattering of HCN emission from the core has prob- 
ably lead to underestimation of the excitation temperature 
and column densities derived using the LTE method in 
Sect. 3.2. For example, an increase of 1 K in Tex, would 
mean 30 % larger column densities, and 3 K would corre- 
spond to a 100% increase in the column density. Quantita- 
tive estimates of the scattering effects are, however, beyond 
the scope of this paper, and here we merely note this addi- 
tional source of error. 

5 Summary 

A mapping in the 7=1-0 rotational lines of HCN, 
HNC and HN'^C shows a compact structure in the di- 
rection of the protostellar core Cha-MMS 1 . The HN'^C 
map agrees best with the previous 1.3 mm dust contin- 
uum map of jR eipurth et al. 1996| , and the HC3N map of 
fKontinen et al. 20001 . 

The column density determinations of H'^CN, HN'^C, 
H'^NC and NH3, and the results of the quoted two stud- 
ies allow us to estimate that the fractional HCN and HNC 
abundances are 7 10"'" and 2 10"^, respectively. These 
are slightly larger than that of NiH^, but clearly smaller 
than the NH3 and HC3N abundances (3 10 and 1 10 ^ 
respectively). In accordance with recent results in other 
dense cores we find that NH3 is about 100 times more 
abundant than N2H^. The kinetic temperature derived 
from the NH3 lines is about 12 K. The derived frac- 
tional abundances conform with the previous suggestion of 
llKontinen et al. 20001 . that Cha-MMS 1 represents an ad- 
vanced chemical state, although it has a large HC3N abun- 
dance, which usually is thought to indicate early times. 
This situation is likely to have developed under the circum- 
stances of CO depletion f llRuffle et al. 19971 V 

The following column density ratios are found: 
HN'^C/H'^CN ~ 3 - 4 and HN'^C/H'^NC ~ 7. Using 
the latter ratio, and the assumption based on the modelling 
results of FT erzieva & Herbst 200 01. we estimate that the 
[HNC]/[HN'^C] ratio is between 30 and 40. As there is 
no HC'^N column density determination available towards 
this source, the possibility that differential '^C fractionation 
causes the observed relatively high [HN'-'C]/[H'-'CN] ratio 
cannot be entirely ruled out. 

The HCN and HNC lines are not much brighter than 
those of the '^C substituted species, but still show no 
absorption dips. These lines are probably weakened 
by scattering by HCN and HNC molecules residing in 
less dense gas around the core. As suggested by 
BCernicharo et al. 19841 . the same process can cause the 
so called HCN hyperfine anomalies in dark clouds, and is 
probably also responsible for the 'too' bright F = - 1 
component of HCN(7 =1-0) towards the centre of Cha- 
MMS 1 . The evidence for scattering and the relatively low 
HCN and HNC column densities derived towards the core 
suggest that the HCN and HNC distributions are less cen- 
trally peaked than those of NH3 and N2H+, which would be 



in agreement with the model of IIAikawa et al. 20051 for a 
rapidly collapsing core. 
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